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Abstract 

In ballistic deposition (BD), (d+ l)-dimensional particles fall sequentially 
at random towards an initially flat, large but bounded d-dimensional surface, 
and each particle sticks to the first point of contact. For both lattice and 
continuum BD, a law of large numbers in the thermodynamic limit estab- 
lishes convergence of the mean height and surface width of the interface to 
constants h(t) and w(t), respectively, depending on time t. We show that 
h{t) is asymptotically linear in t, while w{t) grows at least logarithmically in 
t when d = 1. We also give duality results saying that the height above the 
origin for deposition onto an initially flat surface is equidistributed with the 
maximum height for deposition onto a surface growing from a single site. 



1 Introduction 

Scientific interest in growth processes associated with the deposition of particles on 
surfaces is considerable; see Barabasi and Stanley [3], Cumberland and Crawford [5], 
Vicsek [21]. One family of deposition models involves (d + l)-dimensional particles 
which rain down sequentially at random onto a <i-dimensional substrate (surface); 
when a particle arrives on the existing agglomeration of deposited particles, it sticks 
to the first particle it contacts, which may result in lateral growth and 'overhangs' (if 
it does not contact any previous particle, it sticks to the substrate). This is known 
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as ballistic deposition (BD), and one reason for studying it is as a more tractable 
modification of diffusion limited aggregation (see Atar et al. [1]). 

The physical sciences literature concerned with ballistic deposition in both lattice 
and continuum settings is extensive (see [3] for an overview). As well as numerous 
simulation studies dating back to Void [22], this literature contains analysis by means 
of scaling theory (Family [6], Family and Vicsek [7], Kardar et al. [10]). However, 
these arguments have not been made rigorous from the mathematician's point of 
view; see, e.g., page 56 of [3]. The rigorous mathematical literature is much less 
extensive, but see [1, 16, 17, 19]. The present article builds on the approach of 
Penrose and Yukich [16, 17]. 

In all versions of BD considered here, the substrate is the set W 1 x {0}, identified 
with M. d , or is some sub-region thereof (denoted Q). Thus the substrate is assumed 
initially flat. All particles are incompressible (d + l)-dimensional solids (typically 
balls or cubes). Particles arrive sequentially at random positions in R d . When a 
particle arrives at a position x G M d , its centroid (or some other specified point 
in the particle identifying its location) slides instantaneously down the ray {x} x 
[0, oo), starting from infinity, until the particle hits a position adjacent to either the 
substrate or a previously deposited particle, at which point its motion stops and it 
is permanently fixed. The difference between lattice and continuum models is that 
in the lattice model the positions at which particles arrive are restricted to be in the 
integer lattice Z d (embedded in R d ). 

One is interested in the height and width (roughness) of the interface consisting 
of 'exposed' particles that are 'visible from above'. Loosely speaking, the height and 
width are interpreted here as the sample mean and standard deviation of the heights 
of exposed particles. Let W tjTl denote the width for deposition onto a d-dimensional 
box of volume n, running for "time" (i.e., average number of particles deposited 
per unit volume of the box) t. Scaling theory [6, 7, 10] predicts, and subsequent 
experimental and theoretical studies seem to confirm (for an overview, see [3], [21]) 
that there exist a roughness exponent a and a growth exponent f3, such that the 
surface width is governed by the dynamic scaling relation 

W ttn » n a f(t/n a ^) (1.1) 

where the scaling function / satisfies f(x) oc x 13 , for x <C 1, and f(x) « C for x ^> 1. 
If the scaling theory is correct, it implies that 

W tiH oct 13 t < n; W tyn ocn a , t > n. 

Physicists believe (see e.g. p. 56 of [3]), that dimension d = 1 one has the 'exact' 
values a — 1/2 and /3 — 1/3 (higher dimensional values of a and (3 are not known). 
Validating the above theory rigorously is a challenge for mathematicians; in the 
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present work we take some steps in the direction of rigorously analyzing the regime 
with t Cn. 

In this paper we consider both lattice and continuum BD models. As in [17] we 
first take a thermodynamic limit (deposition onto an infinite surface). In this limit 
we obtain expressions for the limiting height and width, as a function of 'time' (the 
mean number of particles deposited per unit area). The next step is to examine the 
growth of these functions with time, and in doing so we go beyond the continuum 
analysis in [17]. We show that the limiting height grows asymptotically linearly with 
time and that in one dimension, the limiting width grows at least logarithmically. 

The limits are taken in the opposite order in [1] , while the time-parameter and the 
dimensions of the surface are simultaneously re-scaled in Seppalainen [19]. Gravner 
et al. [8, 9] provide detailed results on a one dimensional discrete-time growth model 
sharing some features with BD, which are consistent with the belief that (3 = 1/3. 

2 Lattice ballistic deposition 

We first consider a class of lattice ballistic deposition models, in which all particles 
are assumed identical. Let denote the origin in Z d . Specify a displacement function 
D : 7h d i— > [—00,00) with the properties that (i) -D(O) = 1, and (ii) the set Af := 
{x G Z d : D(x) 7^ —00} is finite but has at least two elements (one of which is the 
origin). For x G Z d let Af x := {x + y : y G Af} and let Af* := {x - y : y G A/"}. The 
set Af is a 'neighbourhood' of the origin and Af x is the corresponding neighbourhood 
of x. The idea of a displacement function is that if a particle arrives at y G Af x , then 
it cannot slide down the ray y x [0, 00) below the position at height D(y — x). 

The substrate is represented by a finite subset Q of 7h d with \Q\ elements. At 
each site x G Q, particles arrive at times forming a homogeneous Poisson process of 
unit rate, independently of other sites. We consider two alternative measures of the 
height of the interface at site x, the last-arrival height Ct,Q( x ) an d the next-arrival 
height r) tj Q(x). The latter is defined in terms of the former by 

Vw (x) := max{&, Q (y) + D(x - y) : y G Z d }, (2.2) 

where we say — oo + rr := —00 for x G R, so in fact r) t Q(x) = max{£ t Q(y)-\-D(x — y) : 

//■■• A';), 

The evolution of 6,q(') proceeds as follows. Assume £o,q(z) = for all z E Z d , 
and as a function of t, assume £,t,Q(z) is right-continuous and piecewise constant 
with jumps only at the arrival times of the Poisson process of arrivals at site z. If a 
particle arrives at site z at time t, then the (last-arrival) height at site z is updated 
to the next-arrival height immediately before time t, i.e. we set 

&,q(x) = Vt-,Q( X ) 1 ™Vs,q{x), 
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while the last-arrival heights at other sites remains unchanged, i.e., Ct,Q(y) — £t-,Q(y) 
for y ^ z. 

Special cases include the so-called nearest-neighbour (NN) and next-nearest 
neighbour (NNN) models [3]. In the NN model, one takes M = {z E Z d : \\z\\i} < 1} 
(i.e., the origin together with its lattice neighbours), and the displacement function 
D is given hy D(x) = for x E Af \ 0; this is the version of ballistic deposition 
considered in [19]. In the NNN model, one takes Af = {z E Z d : H^Hoo < 1} (i.e., 
diagonal neighbours are included) and takes D(x) = 1 for all x E Af; this is the 
version considered in [1]. 

Define the mean height functional £ t g and width functional W t> Q to be the sample 
mean and sample variance, respectively, of the heights at time t, i.e., set 

zeQ y xeQ 

The mean height £ t q is a measure of the average amount of empty space under the 
surface, while the surface width functional W t> Q is a measure of the roughness of the 
interface. 

We consider a particular limiting regime. First we take Q to be large with t fixed 
(deposition for a finite time onto a very large surface), and then we would like to 
take the large-time limit. 

The large-Q limit of £ t Q and W tj Q is best described in terms of deposition onto 
an 'infinite substrate' represented by the whole of Z d . Let £t(z) De the height above 
site z of this infinite interface at time t > 0. Assume again that £t(z) — for all 
z E Z d , but now assume particles arrive as independent Poisson processes at all 
sites in Z d . With this as the only difference from the description of £*,q(')> let the 
updating rules for £ t (-) be just the same as before. Also, define the next-arrival 
height r)t(x) in an analogous manner to (2.2). For the infinite interface process we 
need to check that no 'explosions' occur; our first result does this and more. 

Proposition 2.1 For all t E (0, oo), the values of £ t (0) and r) t (0) are almost surely 
finite, and for all k E N, it is the case that E[(£ t (0)) k ] = 0{t k ) and E[(r] t (0)) k ] = 
0{t k ) as t -> oo. 

We can interpret 6(0) as the height of a 'typical' point in the infinite interface. 
The next result is a thermodynamic limit and shows that in the large n limit, the 
height functional and squared width functional converge to the mean and variance, 
respectively, of the 'typical' height 6(0) (analogously to results in [17] for continuum 
BD). 
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Let (Q n ,n > 1) be a sequence of finite subsets of Z d . Let dQ n denote the set of 
boundary sites in Q n , and set diam(Qri) := max xye Q n \x — y\. Assume that 



G Q n for all n; (2.4) 

liminf(Q„) = Z d , i.e. U~ =1 n™ =n Q m = Z d ; (2.5) 

\dQ n \/\Q n \ -> as n -> oo, (2.6) 

|<5n+i| > |Qn| for all n; (2.7) 

diam((5„)/|(5n| — * as n — > oo. (2.8) 

For example, Q n could be a lattice box of side n centred at the origin. 
Proposition 2.2 For allt G (0, oo), 

lt, Qn £ &(0)] := M*) a « " - ^ (2-9) 

^ t 2 Qn Var[6(0)] := «; 2 (t) as n - oo. (2.10) 



It is likely that with the Poisson arrivals processes at all sites in Z d all put on the same 
probability space, (2.9) and (2.10) also hold with almost sure (a.s.) convergence; 
later we shall give analogous results for continuum BD with a.s. convergence. 

It is possible to give central limit theorems associated with the above laws of 
large numbers. For continuum BD, such results have been given in [17]. Similar 
arguments apply in the lattice case under consideration here, where one can use the 
general lattice central limit theorem of Penrose ([15], Theorem 3.1). 

It is of great interest to estimate the limiting constants h(t) and w(t) in Propo- 
sition 2.2, and especially, to understand the growth of h(t) and w(t) as t becomes 
large. Our main results are concerned with this. Heuristically, one expects h(t) to 
grow linearly in t since the expected height should vary directly with the deposition 
intensity t. The next result demonstrates this, and more. 

Theorem 2.1 There is a constant p\ e (0, oo) such that t~ l h(t) — > p\ as t — > oo. 

Moreover, for any p G [1, oo) we have the LP convergence 

£ _1 &(0) -> pi as t -> oo. (2.11) 

In the special case of the NN model, it can be deduced from Theorem 1 of 
Seppalainen [19] that (2.11) holds with almost sure convergence. Our approach is 
somewhat different from that of [19], and is needed for subsequent results. 

As for the width, the scaling theory mentioned in Section 1 predicts that w 2 (t) = 
0(t 2/3 ). Even without scaling theory, one expects at least that w 2 (t) = 0(t) on 
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the basis of simulations of Zabolitzky and Stauffer ([23], p. 1529), and also on 
the following heuristic grounds. If Af = {0} then the heights £, tin (x),x G Q n are 
independent Poisson variables so that w 2 (t) = t. If J\f ^ {0} so as to give non-trivial 
interactions, these interactions should have a 'smoothing' effect so that w 2 (t) should 
not be any bigger than in the case Af = {0}. 

Rigorous analysis of the large-t behaviour of w(t) appears to be difficult: the 
following result makes a start. 

Theorem 2.2 For any d, it is the case that 

liminfw 2 (t) > 0, (2.12) 

t^oo 

and if d = 1 and Af is a lattice interval, then 

liminf(w 2 (t)/logt) > 0. (2.13) 

t— >oo 

The proof of Theorem 2.2 uses a duality relation between the ballistic deposition 
process and a dual BD process, denoted (£, t (x),x G 1< d )t>o with corresponding next- 
arrival process (fjt(x),x G Z d ) t >o, defined in an identical manner to the original 
BD process £t{x) and next-arrival process r) t (x), except that now one uses the dual 
displacement function D given by D(x) = D(—x), x G Z d , and takes as initial 
configuration a single particle at height at the origin, i.e., one takes 

i (x) = i if x = ° 

^ — oo otherwise. 
We also define a further process (fj t (x),x G Z d ) t > , by 

fjt(x) = rj( t -T)+(x) 

where T is exponentially distributed with mean 1, independent of the process 
(rj s (x),x G Z d ) s > . In other words, the process (f) t (-),t > 0) is obtained by waiting 
an exponentially distributed amount of time before 'kicking off' the dual next-arrival 
process f). We shall refer to f] as the delayed dual BD process. 

Theorem 2.3 Let t > 0. Then the distribution of i] t (0) is the same as that of the 
maximum depth sup zeZ d fj t {z) in the dual BD process. The distribution o/^(0) is 
the same as that of the maximum depth sup^ gZd fj t (z) in the delayed dual BD process. 

Theorem 2.3 shows that w 2 (t) = Var(sup^ eZ d fj t (z)). The dual BD process is a 
random interface growing from a single seed; in this it resembles the classical first 
passage percolation model. Our proof of (2.13) uses ideas from the proof by Pemantle 
and Peres [12] of an analogous logarithmic lower bound for the variance of first 
passage times. Further progress in estimating the growth rate of the variance for 
first passage times has proved elusive (see e.g. [2, 11]); by analogy, the same could 
well be true in the case of w 2 (t). 
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3 Continuum ballistic deposition 



We consider a continuum ballistic deposition model, defined as follows. The sub- 
strate Q is a Borel-measurable region of R d (for example, a cube of side n) and \Q\ 
denotes the Lebesgue measure of Q. Particles are assumed to be (d+ l)-dimensional 
Euclidean balls of possibly random independent identically distributed radii which 
are uniformly bounded by a finite constant -R max (in fact the results presented here 
could easily accommodate a more general class of random shapes, such as convex 
shapes of uniformly bounded diameter). Let F denote the common distribution 
function of the radii of incoming particles; assume that F(0) = and F(R max ) = 1. 

Each incoming particle arrives perpendicularly to the substrate Q x {0} and 
sticks to the first previous particle it encounters, or to the substrate if it does not 
encounter any previous particle. In other words, its motion stops when it encounters 
a previous particle or the substrate, and remains stationary thereafter. 

For simplicity we assume Q is given by one of the sets Q n {n > 1) defined by 

Qn ■= Qn © [0, l) d = {x + y : x G Q n , y G [0, l) d }, 

with the sequence of finite sets Q n C Z d assumed to satisfy the conditions (2.4) - 
(2.8). 

Let V be a homogeneous Poisson point process of unit intensity in M. d x [0, oo), 
each point carrying a mark with distribution F. Denote by Vt, n the restriction of V 
to Q n x [0,t], and denote by V t the restriction of V to R d x [0,t\. 

Represent points in V by (X, T), where X e R d denotes the spatial location 
(center) of the incoming particle and T its time of arrival. Given n, t, the BD process 
driven by Vt, n is defined as follows. The spatial locations of incoming particles are 
given by the spatial locations of the points of Vt, n , the order in which they arrive 
is determined by the time-coordinates of these points (i.e., they arrive in order of 
increasing time-coordinate), and their radii are given by the marks the points carry. 

Let A t) n denote the agglomeration of particles for the BD process driven by 7\ n , 
together with the substrate Q n x {0} (a subset of R d+1 ). For each x G Q n let H t ^ n (x) 
denote the height of the interface above x, i.e., let 

H tin (x) := sup{/i : (x, h) G A t>n }. 

The bold line in Figure 1 represents the graph of the function if ti „(-). 

We define the average height H t ; n and width V^ n of the interface at time t as 
follows. We set H t , n to be the mean of the function H tjn (x),x G Q n , and V tjn to be 
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Figure 1: The mean value of the function whose graph is given by the bold arcs is 
H t n , and its variance is W^ n . The horizontal line represents Q n x {0} 



the root-mean-square deviation of this function from H t n . That is, we set 



H t , n -=\Qn\ 1 / H ttn (x)dx] 



V t , n ■= \l\Qn\~ 1 I (H t>n (x) - H t , n )*dx. 

Qn 



(3.1) 
(3.2) 



These definitions of height and width are slightly different from those used in Penrose 
and Yukich [17] but no less natural. 

The following result gives meaning to the height of the interface above an infinite 
substrate. 



Proposition 3.1 For allt > 0, x e W 1 the limit 

H t (x) = lim H t , n (x) (3.3) 

n— >oo 

exists almost surely and is almost surely finite. Also, the distribution of H t (x) does 
not depend on x. Moreover, for all k G N it is the case that E[(H t (0)) k ] = 0(t k ) as 
t — > oo. 



Thus the infinite-substrate height function is given by H t (x),x e M. d . The fol- 
lowing thermodynamic limits are variants of results in [17]. They are continuum 
analogues to Proposition 2.2. In particular, they show that the variance of the ran- 
dom variable H t (0) (height above a typical point of the infinite substrate) is the 
large-n limit of the width functionals V t 2 n (sample variance of the heights above a 
bounded substrate). 
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Proposition 3.2 For allt G (0, oo), we have the almost sure convergence 

\im(H t>n )=E[H t (0)]; (3.4) 
limlf n = Var[tf t (0)]. (3.5) 

n— »oo 

Now we give a continuum analogue to Theorem 2.1 

Theorem 3.1 There is a constant p 2 G (0, oo) such that t^ElH^O)] — > p 2 as 
£ — > oo . Moreover, for any p > 1 , 

f -1 i?t(0) -> p 2 in L p as t -> oo. (3.6) 

The next result is the continuum analogue to Theorem. 2.2. 

Theorem 3.2 For any d 7 i£ is the case that 

liminf Var[# t (0)] > 0, (3.7) 

and if d = 1, then 

liminf(Var[# t (0)]/log£) > 0. (3.8) 

Also of independent interest is the duality result given by Proposition 5.1 below, 
which is the continuum analogue of Theorem 2.3. 

4 Proofs for lattice BD 

For x G Z d , let S\(x), S^x), S^x), . . . denote the ordered arrival times of the Poisson 
process at site x. The proofs for lattice BD are based on a directed graph repre- 
sentation. For each t > 0, define a directed graph Q t and a directed graph Q t , both 
with vertex set V t defined by 

V t := {(z, Si(z)) : z G Z d , i G N, S^z) < t} U {(x, 0), (x, t); x G Z d }. 

Informally, each Poisson arrival at site x G Z d before time t is represented by a point 
(x, T) G Z d x [0, oo). Two points (x, T) and (y,U) of Vt are joined by a directed 
edge in Q t from (x, T) to (y, [/) if T < U and y G A/" x . They are joined by a directed 
edge in Q t from (x, T) to (y, U) if T < U and y G M*. 

A pa£/i in (/ t is a sequence 7r of vertices in V t denoted (x ,T ), . . . , (x n ,T n ) say, 
such that T = and T n = t and for i = 1, 2, . . . , n there is an edge of Q t from 
to We say the path starts at (x ,0) and ends at (x„,£). The 

length of the path is n. The height of the path, denoted /i(7r), is X^=i -^(^j — ^j-i)- 
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A path in Q t is denned similarly, except that now the requirement is that for 
each i there is an edge of Q t from to (xj,Tj), and the (dual) height h{ii) 

of a path 7i = ((xi, T;))™ =0 in Q t is given by £)™ =1 D(xj - Xj-i). 

The skeleton of a path 7r in Q t is its projection onto Z d , i.e., the sequence 
(x , • • • , x n ). Given a sequence (x ,...,x n ) G (Z d ) n+1 with x^ — Xj_i G A/" for 
1 < i < n, a maximal path with skeleton (xi, . . . , x n ) is a path ((x , T ), . . . , (x n , T n )) 
in the graph Q t , with T = and T n = t, and with the property that for 1 < i < n, 
there are no Poisson arrivals at x n in the time-interval (Tj_!,Tj). 

Lemma 4.1 Let £ > and z G Z d . Suppose the maximum length of all paths in V< 
ending at (z,t) is finite. Then r) t -(z) := lim s f t ?7 s (z) is g'wen &y 

r)t-(z) = sup{/i(7r) : 7r a path in C/ t ending at (z,t)}, (4.1) 

and f) t -(z) := lim s -f t 77 s (z) gwen &y 

fjt-(z) = sup{h(n) : 7r a path in Q t starting at (0,0) and ending at (z,t)}, (4.2) 

with the convention that the maximum of the empty set is — oo. 

Proof. If 7i is a finite path ending at (z,t), then we assert that r) t -(z) > h{ji). 
Indeed, if tt — (xj,Tj)™ =0 then by monotonicity of the processes £ t and rj t , for each i 
with 1 < i < n we have 

fri(Xi) > VTi-(Xi) > fri-Axi-l) + D ( x i ~ x i-l) 

so that 

Vt-(z) > ^T n ^(Xn-l) + D(x n - X n _i) 

n 

> D(xi - Xi_i) = h(n). 
i=i 

Conversely, there is at least one path of height at least rjt-(z) that ends at (z,t). 
This is proved by induction on the maximum length of paths ending at (z,t); it 
is clearly true when this maximum path length is 1; suppose it is true when the 
maximum path length is in the range {1,2, ... ,k}. Now suppose the maximum 
path length is k + 1. Then 

r] t -(z) = max{&_(y) + D(z - y) : y G M* z } 

so that for some y* G N* (i.e., with z G N y *) we obtain 

ifc-(«)=6(y*)+I>(«-y*), 
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and if L denotes the last Poisson arrival time before time t at site y*, we have 



rk-(z)=Z L (y*) + D(z-y*) 
= VL _(y*)+D(z-y*). 

Each path in Q L ending at (y*, L) has length at most k, since otherwise there would 
be a path through (y*, L) to (z, t) of length greater than k + 1. Hence by the induc- 
tive hypothesis, there is a path in Ql ending at (y*,L), of height at least i]L_(y*), 
and hence by appending (z, t) to this sequence one obtains a path in Q t , ending at 
(z, t), of height at least r] L ^(y*) + D(z — y*). This completes the induction and hence 
the proof (4.1). The proof of (4.2) is similar, g 

Proof of Proposition 2.1. Since ^(0) < i] t (0), it suffices to prove that i] t (0) is a.s. 
finite and E[r] t (0) k ] = 0{t k ) as t -> oo. Note that 7ft (0) = rjt-(0) a.s. Also, let 
-Dmax : = m ax zeZ d D(z) denote the maximum value of the displacement function. By 
our assumptions on this function, we have 1 < -D max < oo. 

By Lemma 4.1, if ?7t-(0) > m then there must be a path in Q t of height at least 
m, and hence of length at least m/D ma , x , that ends at (0,t). Hence, if ?7t-(0) > m 
then there is a path in Q t of length at least m/D mSuX that is maximal for its skeleton 
and ends at (0, t). 

For any given sequence (xo, x±, . . . , x n ) with Xi — G H for i = 1, . . . , n, the 
probability that there exists a maximal path in Q t with skeleton ( ) and 

with all arrival times less than t, equals than the probability that the sum of n — 1 
independent exp(l) variables is less than t, which is the same as the probability that 
Po(£) > n — 1, where Po(t) denotes a Poisson variable with mean t. Hence, for any 
y >te 2 + 1, by e.g. Lemma 1.2 of [13], 

P[Vt(0) > yD max ] = P[vt-(0) > yD max ] < \Af\ y P[Po(t) >y-l] 

< |Ar|^exp(-((y-l)/2)log((y-l)/t)) (4.3) 



which tends to zero as y — > oo. Hence r] t (0) is almost surely finite. Also, for k E N, 





V^(o)y" 


POO 

= P 




\ -^max / 


Jo 



max 



Set c := (2|A/"|) 2 , split the region of integration into w < (ct < l) k and w > (ct+l) k , 
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and change variables to y — w 1 ^ in the second integral to obtain the estimate 



and the last integral is finite, so that t k E[r) t (0) k ] is bounded, g 

Proof of Proposition 2.2. The idea is to apply Theorem 3.1 of Penrose [15]. Let B 
denote the collection of all non-empty finite subsets of Z d . Let X = (X x , x G Z d ) be 
a family of independent homogeneous Poisson processes of unit intensity. For Q G B, 
assume the evolution of ^ t ,Q is governed by the Poisson processes (X x , x G Q). Then 
(£t,Q(x),Q £ B, x G Q) is stationary £>-indexed summand in the sense of section 3 
of [15]. Also, the proof of Proposition 2.1 shows that 



For any B- valued sequence (B n ,n > 1) with liminf(_B n ) = Z d , we have as n — > oo 
that £t,B n (0) — > 6(0) almost surely (see the proof of Proposition 2.1, or Lemma 5.1 
of [15]), and hence in L? by (4.4). Hence, for any such sequence (B n , n > 1) we have 
the L 1 convergence 



By (4.5), the first part of Theorem 3.1 of [15] is applicable to the stationary B- 
indexed summand (Ct,Q( x ),Q G £>, x G Q), so that eqn (3.3) of [15] gives us (2.9). 

Moreover, since (4.5) holds with L? convergence as well, by changing L 1 estimates 
to L? estimates throughout the proof of eqn (3.3) of [15], we may deduce that result 
in this case with L 2 convergence, i.e. (2.9) holds with L 2 convergence as well (the 
proof in [15] in its turn uses a multiparameter L 1 ergodic theorem quoted from [14], 
but this is also easy to extend to L 2 convergence in the present setting, using (4.4).) 

To prove (2.10), first expand the sum of squares in (2.3) to obtain 




< (ct + l) k + ke 1 / 2 




sup{El£ t< Q(x) 4 ] : Q G B,x G Q} < oo. 



(4.4) 



6a(0)-6(0); 
6a(0) 2 ->6(0) 2 . 



(4.5) 
(4.6) 




(4.7) 
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By (4.6), the first part of Theorem 3.1 of [15] is applicable to the stationary B- 
indexed summand (Ct,Q( x ) 2 , Q G B, x G Q), and this shows that the first term in the 
right hand side of (4.7) converges in L 1 to E[^ t {0) 2 ]. 

Since (2.9) holds with L 2 convergence, it follows that the second term in the 
right hand side of (4.7) converges in L 1 to (E[£ t (0)]) 2 . Combining these limiting 
results in (4.7), we obtain (2.10). | 

Proof of Theorem 2.3. The idea of the proof is a form of time-reversal of the 
graphical representation. Let ip t : Z, d x [0,t] — > Z d x [0,t] be defined by 

ipt((z,s)) = (z,t- s). 

Let V t := -0 t (Vt) , and observe that the point set V t has the same distribution as V t 
by the properties of Poisson point processes. 

Let Qt be defined in the same manner as Qt but on the vertex set Vt instead of 
Vt- Then each path in Q t ending at (0, t) corresponds to a path with the same height 
starting at (0, 0) in Q t ; the correspondence is obtained by reversing the sequence of 
vertices and then applying the mapping ip t to each vertex in the sequence. 

Note that r] t (0) = i] t -(0) with probability 1. By (4.1), i] t (0) is the greatest height 
of all paths in Q t which end at (0,t). By the correspondence described above, this 
is precisely the same as the greatest height of all paths in the graph Q t which start 
at (0,0). 

Since the point processes Vt Vt have the same distribution, it follows that i] t (0) 
has the same distribution as the greatest height of all paths in the graph Q t which 
start at (0, 0). Hence by (4.2), r) t (0) has the same distribution as sup zgZ d fjt-(z) and 
hence the same distribution as sup 2gZ d fj t (z). 

Let L be the last arrival time at before time t (or L — if there are no arrivals 
at before time t). Then t — L has the distribution of min(T, t) where T is expo- 
nential with mean 1, and hence L has the distribution of (t — T) + . Given L with 
L > 0, the distribution of ^_(0) is the same as that of sup zgZ d f)L{z), by the same 
argument as above and the Markov property of the time-reversed Poisson process. 
Whenever L > we have &(0) = ^l-(O), and if L — then &(0) = almost 
surely. Combining these observations shows that 6(0) has the same distribution as 
sup 2gZ d fj(t-T)+(z), as asserted, g 

From now onwards, we shall assume the delayed dual BD (next-arrival) process 
fjt(z) and also the associated last-arrival process £t( z ) is defined in terms of the 
Poisson arrival times (Si(x),i > l,x G Z d ), as follows. For t < Si(0) we put 
fjt(z) = and £ t (z) = — oo for all z G Z d . Then we put 

€ (z) ■= { if z = 

x ^ ' ^ — oo otherwise. 
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and define f]Si(o) in terms of ^(o) in the usual manner as given at (2.2). Then we 
allow the evolution of (Ct,Vt)t>Si(o) to follow the usual rules driven by the Poisson 
arrivals {Si(z) : Si(z) > Si(0). Then fj t (x),x > constructed in this manner clearly 
follows the desired evolution prescribed in Section 2, with Si(0) used as the initial 
'kicking off time'. 

For t > 0, u > 0, let D t be the depth (i.e., the maximum next-arrival height) 
of the delayed dual BD process fj t at time t, and let T{u) be first passage time to 
depth u for the delayed dual BD process, i.e. let 

D t := sup{fj t (z)}- T{u) := m£{t : D t > u}. 
Lemma 4.2 There is a constant p e (0, oo) such that 

lim = p, a.s. (4.8) 

and 

lim — = p~\ a.s. (4.9) 

t^co t 

Proof. First we verify that E[T(1) 2 ] is finite. Note that fj t (0) > N t — 1, where here 
N t denotes the number of arrivals at up to time 1. Hence, 

P[T{1) >t}< P[fj t (0) < 1] < P[N t < 2] = P[X + Y > t] 

where here X and Y denote independent exponential random variables with unit 
mean, representing the first two inter-arrival times at 0. Hence, 

E[T(1) 2 } < E[(X + Y) 2 } < oo. 

Next, we assert that T{u) is distributionally subconvolutive, i.e. for u, v > we 
have 

F T (u+v) > F T(U) * F T(V) . (4.10) 

To see this, let X{u) G Z d be chosen (in an arbitrary way if there is more than one 
choice) so that fiT( u ){X{u)) > u (by definition such X(u) exists). Let T(u) + T* 
be the time of next Poisson arrival after T(u) at site X(u), and let {is,Vs)s>o be a 
version of the BD process with displacement function D, with initial profile 




iix = X(u) 
— oo otherwise 
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and driven by Poisson arrival times {S*(x)} given for each x E Z d by 



{S*(x)} = {Si(x) - T{u) - T* : Si(x) > T{u) + T*}, 

where {Si(x)} are the arrival times driving the original dual BD process (£t,Vt)- 
Let T**(v) be the first time the process (rj*) s > achieves a depth of at least v, i.e. 

T**{v) = inf{t > : sup(fj*(x)) > v}. 

Then T* + T**(v) has the same distribution as T(v), and is independent of T{u). 
Also, since £t(u)+t*(X(u)) > u, the depth at time T(u) + T* + T**(v) is at least 
■u + f , i.e. 

SUp(r/T( u ) +T *+T"(,;)(^)) > « + V, 

so that T(u) +T* + T**(v) > T(u + u). Combining these facts gives us (4.10). Since 
the variables (T(u),u > 0) are also monotonically increasing in u, we can apply the 
Kesten-Hammersley theorem ([20], page 20) to obtain (4.8), with < p < oo. Also, 

AAA A + 1 /x 

> _L > L_ x - — - 4.11 

T(A) " * " A + 1 T(A + 1) 

and since A — *■ 00 almost surely this with (4.8) yields (4.9) provided p > 0. 

If p = then by (4.8) and the second inequality of (4.11) we would have D t /t — > 
oo almost surely, so that E[D t /t] — > oo, and so by Theorem 2.3, E[^ t (0)/t] — > oo. 
This contradicts Proposition 2.1, and hence p > as asserted, g 

Proof of Theorem 2.1. Set pi := p _1 , with p as given in Lemma 4.2. By Theorem 
2.3, A has the same distribution as r) t (0), and so by Proposition 2.1, for any p > 1 
the pth moment of (D t /t) is bounded uniformly in t. By (4.9), (A A) converges 
almost surely to pi, and by the moment bound the almost sure convergence extends 
to convergence in pth moment for any p > 1. Since A has the same distribution as 
r)t(0), this convergence in pth moment also holds for 77* (0) . | 

Proof of Theorem 2.2. To prove (2.12) consider the event that (i) there are no 
arrivals in M\ {0} between times t — 1 and t, and (ii) there is at least one arrival 
at in the time-interval (t — l,t]. This event has the same non-zero probability 
for all t > 1. Conditioned on this event, and on everything before time t — 1, the 
conditional variance of £t(0) is the variance of the number of Poisson arrivals at 
in the time-interval (t — l,t], i.e. the variance of a Poisson variable with unit mean 
conditioned to take a value of at least 1. This variance is a strictly positive constant 
and (2.12) follows. 
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Now suppose d — 1 and ./V is a lattice interval. To prove (2.13), we consider the 
delayed dual BD process {£,t,Vt)t>o- 111 this process, we denote by accepted arrival 
an arrival time T = Si(z) such that £t(z) > —oo. Enumerate the set of all accepted 
arrival times (for all sites) in increasing order as Ti, r 2 , . . .. Let N t be the number of 
accepted arrivals up to time t, i.e. 

N t := sup{n : r n < t}. 

Let I t be the size of the interface at time t, i.e., the number of sites z E Z with 
fjt(z) > — oo. For n > 1, let Y n := I Tn be the size of the interface after n accepted 
arrivals, and set Yq := 1. Since (N t ,t > 0) is a Poisson counting process with its 
'clock' running at speed I t , we have that 

lim *_ = 1, a.s. (4.12) 

Next, we assert that there is a constant 7 > such that 

lim (V*) = 7, a.s. (4.13) 

To see this, recall that we are assuming here that M is a lattice interval including 
and at least one other element. It is not hard to see that I t must also be a lattice 
interval, and that both the right and the left endpoint of I t follow renewal reward 
processes, where in both cases the the inter-arrival times of the underlying renewal 
process are exponentially distributed and the rewards are uniformly distributed over 
a lattice interval. The assertion (4.13) follows by the Strong Law of Large Numbers 
for a renewal reward process. 

By (4.12) and (4.13) we obtain 

Nt 1 

hm-i = I, a.s. (4.14) 

t— >oo t A 

Let M(u) be the number of accepted arrivals up to time T(u). By (4.14) and (4.8), 
as u — > 00, it is the case with probability 1 that 

M(u) = N T(U) ~ (i/2)T(u) 2 ~ (ip 2 /2)u 2 . (4.15) 

Let T be the a-algebra generated by the locations in Z of the sequence of accepted 
arrivals. Conditional on JF, the distribution of T(u) is that of the sum of M(u) 
independent exponentials with the jth exponential having mean Y~\. Hence, 

M(u) 

al:=V a r[T(u)\T]=J2 Y A- 
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By definition, N T . = j, and hence by (4.14), j/rf — > 7 /2 so that Tj ~ {2j/^) l l 2 as 
j — > 00, almost surely. Since Yj = I Tj , by (4.13) we obtain as j — > 00 that with 
probability 1, 



Y} ~ 7Tj ~ (27J) 1 / 2 , 
so that by (4.16) and (4.15), as u — > 00 we have that 

<t„ ~ (27) _1 log M(-u) ~ 7 _1 log u, a.s. 



(4.17) 



(4.18) 



The Berry-Esseen theorem, e.g. as given in theorem 5.4 of Petrov [18] or as 
quoted in Chen and Shao [4], says that there is a constant C such that if Xi, . . . Xk 
are independent random variables with mean zero and finite third moments and 
W := ^2i =1 Xi has variance 1, then 



sup{\P[W <x]- S(x)|} < C*]Tp[|X,| 3 ], 



i=i 



where <£> is the standard normal distribution function. 

Let r := 0, and for i > 1 let e$ := Tj — Tj_i. As mentioned earlier, conditional 
on T the are independent exponentials with _E[ej|jF] = Y^. Define 

M(u) 

= ^ J E[|e l - J E[e i |^]|-V] 



Af(tt) 

= (12c- 1 - 2) £ 



(4.19) 



i=i 



since if X is exponential with mean a then E[\X — a\ 3 } - 
Set fi v := E[T(u)\F\. By the Berry-Esseen Theorem, 



sup 



P 



T(u) - n u 



< x\T 



<7„ 



< 



a 3 (12e" 1 -2). 
CO,, 



(4.20) 



By (4.17), the sum Y2i=i ^-1 converges almost surely. Hence by (4.18) and (4.19), 
we can find uq such that for u> uq we have -PL4J < 0.01 where 



C9 U 



> 



0.01 L U K < V / (log«)/(27)} 
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so that A u e T . Then for any y e R, using (4.20) we may deduce that 



= P 



- ^ < y-/^u + o.2( 7 - 1 iog M ) 1 / 2 ^ 



< y + 0.2( 7 - 1 log uy' 2 \AV\ - P[T(u) < y\A c u ] 
T(u) - fj, u y - \i u 



T(u) - n u y - n u 



< P 



< 



+ (0.2)v/2K 



Z — \ n u 



< sup {$(x + (0.2)V2) - $(x)} + 0.02 < 0.22, 

xeR d 

and so for u > u$ we have 

P[T{u) <y + 0.2( 7 - 1 logw) 1/2 ] - P[T{u) <y}< 0.22 + P[A U ] < 1/4. (4.21) 
For t > 0, let z/(t) be the median of the distribution of D t . Then 

P[T(z/(t)) >t] = P[D t < v(t)\ < 1/2 (4.22) 
so that P[T(v(t)) <t}> 1/2, and hence by (4.21), 

P[T(u(t)) < t - 0.2( 7 " 1 \ogu(t)) 1/2 } > 1/4. 
By (4.9), (D t /t) — > in probability, so that (u(t)/t) — > as t — > oo. Hence, 

\ogv(t) ~ logt as £ — > oo. (4.23) 
If T(v(t)) < t - 0.2(7 _1 log z/(t)) 1/2 , then there is a site z*6Z with 

^-0.2(7"! log ^(t)) 1 / 2 !^*) > ^(*)- 

If also at least 0.1 ( 7 _1 logt) 1 / 2 + 1 particles arrive at site z* between times t — 
0.2(7 _1 logi/(*)) 1 / 2 and t, then D t > u(t) + O.l^ 1 logt) 1/2 . By (4.23), the condi- 
tional probability of the second of these events, given the first, tends to 1 as t — > oo, 
so that for large t, 



P[D t > u(t) + O.l^- 1 logt) 1/2 | > 1/8. 



(4.24) 



Moreover, by definition P[D t < v(t)\ > 1/2. Combining this wtih (4.24) yields 
(2-13). B 
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5 Proofs for continuum BD 



In this section, for iGl 1 * and r > we write B r (x) for the closed Euclidean ball of 
radius r centred at x. 

We introduce a dual continuum BD process, denoted H t (x), a defined in an 
identical manner to the original continuum BD process except that now the initial 
profile H (x),x G M. d is given by 

H (x) = { ° if X = ° 

°^ ' — oo otherwise. 

In other words, the interface grows from an initial seed consisting of a single point 
at (0,0). Incoming particles which miss the agglomeration have no effect. 

We adapt the graphical representation to the continuum. Given a locally finite 
marked point set X C M d x [0, oo), we say a sequence of points {X^T^)\ =1 in X 
forms a path in X if T\ < T2 < • ■ • < T^, and for 1 < % < k, the ci-dimensional ball 
centred at Xi (of radius given by the mark of (Xj,Tj)) overlaps the ball centred at 
X i+ i (of radius given by the mark of (X i+ i,T i+ i)). 

Given such a path, we refer to the sequence (Xi, . . . , X k ) as the skeleton of the 
path. Given also x G M. d , we say the path starts near x if x lies in the <i-dimensional 
ball centred at X\ (of radius given by the corresponding mark) and we say the path 
ends near x if x lies in the d- dimensional ball centred at X n . We say the height at x 
of the path is the height of the interface above x for the BD agglomeration onto an 
initially flat surface determined by the finite sequence of incoming (d+l)-dimensional 
balls centred at Xi, . . . , X k (with radii given by the corresponding marks). 

If the path starts near 0, we say the dual height of the path at x is the height of 
the interface above x for the BD agglomeration determined by the finite sequence 
of incoming (d + l)-dimensional balls centred at X±, . . . , Xk (with radii given by the 
corresponding marks), starting initially the initial profile being given by the function 
H (i.e. by a single point at (0,0) rather than an initially flat surface). 

Lemma 5.1 Let t > 0, x G M. d . Then with probability 1: 

(i) There are a.s. only finitely many paths in Vt which end near x; 

(ii) the value of H t (x) is the maximum height at x of all such paths (or is zero 
if no such path exists ); 

(Hi) H t ^ n (x) is the maximum height at x of all paths in Vt which end near x for 
which the skeleton is contained in Q n , (or is zero if no such path exists), and 

(iv) H t (x) is the maximum dual height at x of all paths which start near and 
end near x (or is —00 if no such path exists). 

Proof. Part (i) follows from e.g. corollary 3.1 of [14]. 
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In the continuum BD process, there exists a sequence of particles in the ag- 
glomeration, each particle touching the next one in the sequence, leading from the 
substrate to a particle arriving at time at most t and with (x, H t (x)) on its surface 
in the agglomeration; considering only the particles in this sequence, we have a path 
which ends near x and has height H t (x) at x. Hence the height H t (x) is at most the 
maximum height at x of all such paths. 

On the other hand, inserting extra points into a given (d+ l)-dimensional marked 
point process cannot decrease the height over x of the interface of the corresponding 
BD agglomeration, so for any path in Vt ending near x, H t (x) is at least the height 
at x of the path. This demonstrates part (ii). 

The argument for part (hi) is the same as for part (ii), except that now one 
ignores particles with spatial location outside Q n . 

The argument for part (iv) is the same as for part (ii). g 

Proof of Proposition 3.1. Let x G M. d . By (2.5) and part (i) of Lemma 5.1, there 
exists an almost surely finite random n such that for n > n , every path in V t 
which ends near x has its skeleton contained in Q n . Then by parts (ii) and (ii) of 
Lemma 5.1, we have H t (x) = H tjn (x) for n > n . This demonstrates the first part 
of Proposition 3.1. 

Since Vt is distributionally invariant under spatial translations (i.e., translations 
of R d x R leaving the time-coordinate unchanged), the distribution of H t (x) := 
lim^oo H t>n (x) does not depend on x. This demonstrates the second part of Propo- 
sition 3.1. 

We prove the last part only in the case where R max < 1/2; the more general case 
can then be deduced by some simple scaling arguments which we omit. 

We couple our continuum BD process to a certain NNN lattice BD model, as 
defined in Section 2. Partition M. d into half-open unit cubes, and for iGl' i let Q(x) 
be the the cube in the partition that contains x. Let (£t, f]t)t>o be the coupled NNN 
lattice BD model, in which the arrival times at site z G Z d are given by the time- 
coordinates of the points of V in Q(z) x (0, oo). By the assumption that 2R max < 1, 
it is not hard to see that H t (x) < r] t (z(x)). We can then use Proposition 2.1 to 
deduce that E[(H t (0)) k ] = 0{t k ). ■ 

Proof of Proposition 3.2. By Lemma 5.1 and distributional invariance of Vt under 
spatial translations, there is a function q(t),t > such that q(t) — > as t — > oo and 

P[H t {x) + H t , n (x)] < q(dist(x, dQ n )). (5.1) 

Also, by the argument in the proof of the last part of Proposition 3.1, for any k > 1 
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we have 



sup E[H t , n (x) k ] < oo. 



(5.2) 



n>l,a;€lR d 



By (3.1), (5.2) and Fubini's theorem, 



E[H t , n \ = \Q 



n 



-1 



/ 



n 



EH t ^ n (x)dx. 



(5.3) 



By (5.1) and (5.3), using the condition (2.6) on Q n it is straightforward to deduce 
that EH t n — > EH t (0) as n — »■ oo. 

An alternative representation of H t n is as the sum 



where £((X,T),Vt, n ) is integrated height of the vertically exposed parts (at time t) 
of the ball arriving at (X, T). In other words, taking J(x,t)( x ) to be the indicator 
function of the event that the closure of the ball deposited above location X at time 
T includes a point with with coordinates (x, H ttn (x)), we set 



Using (5.4), we can obtain a law of large numbers with almost sure convergence and 
convergence of the mean, for H tjTl (x) by use of Theorem 3.2 of [16] (this is where we 
use conditions (2.7) and (2.8) on Q n ). As mentioned above, the mean converges to 
EH t (0), so the limiting constant must be EH t (0). This demonstrates (3.4) 

We now prove (3.5). By expanding out the square in (3.2), we obtain the identity 





(5.4) 



(x,T)er t} 





(5.5) 



so that by Fubini's theorem 



E[V t 2 n ] + [EH t , n ] 2 = IQJ- 1 / E[Hl n {x)\. 



Using (5.1), (5.2) and (2.6), it is then straightforward to deduce that 



E[V t 2 n ] + [EH ttn f - E[(H t (0)) 2 ] 



(5-6) 
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By (5.5), we have 

V t 2 n + Hl n = \Qn\~ l £ ^ 2 \(X,T),V t , n ) (5.7) 

{X,T)&Vt,n 

where we set 

^ 2 \(X,T),V t , n )= [ H t , n (x) 2 J {x , T) (x)dx. 

JR d 

Hence we may use Theorem 3.2 of [16] to deduce that the expression (5.7) converges 
to a deterministic limit, almost surely and with corresponding convergence of mean 
(the 'bounded pth. moments condition' in [16] can be shown for this case by argu- 
ments in [16].) Hence by (5.6), the almost sure limit of the expression (5.7) is equal 
to E[(H t (0)) 2 ]. Combining this with (3.4), we obtain (3.5). ■ 

For t > and u > 0, let D t denote the depth (i.e. maximum height) of the dual 
continuum BD model at time t, and let T(u) be the first passage time to depth u 
of the dual continuum BD model, i.e. let 

A := sup (H t (x)); f(u) := M{t > : D t > u}. 

Proposition 5.1 The distribution of H t (0) is the same as that of D t . 

Proof. Fix t > 0, and consider the time- reversed space-time Poisson process with 
time-coordinates transformed by the mapping (X, T) i— > (X,t — T). Under this map- 
ping, any path ((X 1 ,T 1 ), . . . , (X k , T k )) corresponds (by reversing the order of points) 
to a path in the transformed Poisson process, namely ((X k ,t — T k ), . . . , (X 1 ,t — T x )), 
the so-called time-reversed path. If the original path ends near 0, the corresponding 
time-reversed path starts near 0, and the height at of the original path equals the 
dual height over X\ of the time-reversed path, which is the maximal dual height of 
the time-reversed path. 

By Lemma 5.1 (ii), H t (0) > u if and only if there is a path in Vt which ends near 
with height at of at least u by time t, in which case the corresponding time- 
reversed path has maximal height at least u. Hence, H t (0) is the maximal dual 
height of time- reversed paths in Vt starting near 0. Hence by Lemma 5.1 (iv), H t (0) 
is the maximal depth at time t for the continuum ballistic deposition process driven 
by the transformed Poisson process, using initial profile Hq. Since the distributions 
of the original and transformed Poisson processes are identical, H t (0) therefore has 
the same distribution as maximal depth in the BD process generated by the original 
Poisson process with initial profile H . In other words, it has the same distribution 
as Df ■ 

The next result is a continuum analogue to Lemma 4.2. 
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Lemma 5.2 There is a constant p% G (0, oo) such that 



Inn = p 3 , 



a.s. 



and 



lim — - 

t^oo t 



P3 



a.s. 



(5.8) 



(5.9) 



Proof. First we show that T(l) has finite second moment. Choose e\ > such that 
F(e±) < 1/2, and e 2 > such that a ball of radius e\ in contains a rectilinear 
cube of side e 2 , with the same centre. Let Ni(t) be the number of Poisson arrivals in 
Vt having spatial coordinate with norm at most £2/2 and having mark at least 
e\. Then Ni(t) is Poisson with parameter at least e 2 t/2, and Ht(0) > e 2 N 1 (t), so 
that 



P[f{\) >t]< P[H t (0) < 1] < P[Yo{e 2 t/2) < l/e 2 ] = P 



l<i<ri/e 2 l 



where e^e^eg, . . . are independent exponential random variables with mean 2/e 2 , 
representing the inter- arrival times of a Poisson process of rate e 2 /2, and \x] denotes 
the smallest integer not less than x. Hence, 



E[f(l) 2 } < E 




< 00. 



We assert that T(u) is distributionally subconvolutive, i.e. for u, v > we have 



Ff( u+V ) ^ Ff(u) * Ff(v)- 



(5.10) 



To see this, let X(u) G M. d be chosen (in an arbitrary way if there is more than one 
choice) so that Hf,JX(u)) > u (by definition such X(u) exists). Let (H*(x),x G 
>o be a version of the BD process with initial profile 



H*(x) 



if x = X(u) 
—00 otherwise 



and driven by the Poisson process T_f, u JV) fl {R d x (0, 00)), where r_ t denotes the 
shift operator on R d x R mapping each point (x, u) G M. d x R to (x, u — t). 
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Let f*(v) be the first time the process (H*) s > achieves a depth at least v, i.e. 
f*(v) = inf{t > : sup (#;(&)) > v}. 

x£R d 

Then T*(v) has the same distribution as T(v), and is independent of T(u). Also 
the depth at time T(u) + T*(v) is at least u + v, i.e. -Df r^+f*^) > w + t>, so that 
T(-u) + T*(t>) > T(w + v). Combining these facts gives us (5.10). 

Since the variables (T(u),u > 0) satisfy (5.10) and are also monotonically in- 
creasing in u, we can apply the Kesten-Hammersley theorem ([20], page 20) to obtain 
the desired conclusion (5.8), with < p 3 < oo. 

The arguments to show that (5.9) holds and p 3 > are just the same as the 
corresponding arguments in the proof of Lemma 4.2. g 

Proof of Theorem 3.1. The proof is entirely analogous to that of Theorem 2.1, 
now using Propositions 3.1 and 5.1 along with Lemma 5.2. | 

Proof of Theorem 3.2. First we prove (3.7). Let Tt be the a- field generated by 
all arrivals up to time t — 1. Then H t -\(0) is T t - measurable. It is clear that there 
is a constant e > such that 

P[H t (0) = H t ^(0)\T t ] > e a.s. 

and 

P[H t (0) > # t _i(0) + l\T t ] > e a.s. 

and combining these two estimates shows that the conditional variance Vai[H t \Tt\ 
is bounded away from zero. From this we may deduce (3.7). 

The proof of (3.8) is similar to that of (2.13). Now take I t to be the Lebesgue 
measure of the interface, i.e., of the set of sites x E M. d with H t (x) > — oo, and take 
T to be the a-algebra generated by the sequence of locations and marks of accepted 
arrivals in the dual continuum BD process. 

Then it is again the case that conditional on T, the distribution of T(u) is the 
sum of M{u) independent exponentials e±, . . . ,eM( u )] this is because, given T, the 
ej are independent exponentials. Indeed, given the positions of the first j accepted 
arrivals, the distribution of ej is exponential with mean Y~ l ; extra information about 
the location of this arrival and subsequent arrivals does not affect its distribution. 

Using this information, the proof of (3.8) follows that of (2.13) closely, and we 
omit further details. | 
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